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ADDRESS OF SMALL BUSINESS CONCERN: 6751 Rufe Snow Drive. Suite 
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fiscal year, and (2) concerns are affiliates of each other when either, directly or 
indirectly, one concern controls or has the power to control the other, or a third 
party or parties controls or has the power to control both, 

I hereby declare that rights under contract or law have been conveyed to 
and remain with the small business concern identified above with regard to the 
invention described in: 

[X] the specification filed herewith with title as listed above. 
[ ] the application identified above. 
[ ] the patent identified above. 

If the rights held by the above-identified small business concern are not 
exclusive, each individual, concern or organization having rights in the invention 
must fife separate verified statements averring to their status as small entities, 
and no rights to the invention are held by any person, other than the inventor* 
who would not qualify as an Independent inventor under 37 (XF.R. 1.9(c) if that 
person made the invention, or by any concern which would not qualify as a small 
business concern under 37 C.F,R. 1.9(d), or a nonprofit organization under 37 
(XF.R. 1.9(e) 

Each person* concern or organization having any rights In the invention 
is listed below: 

[X] no such person, concern, or organization exists. 

I ] each such person, concern or organization is listed below- 

FULL NAME: 
ADDRESS: 

[ Individual [ ]SmaiI Business Concern [ JNonprofit Organization 
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Separate verified statements are required from each named person, 
Concern or organization having rights to the invention averring to their status as 
small entitle. (37 C.RR. 1 .27). 



I acknowledge the duty to file, in this application or patent, notification of 
any change In status resulting in loss of entitlement to small entity status prior to 
paying, or at the time of paying, the earliest of the issue fee or any marntenance 
fee due after the date on which status as a small entity is no longer appropriate. 
{37 1.28 (b)) 



I hereby declare that all statements made herein of my own knowledge 
are true and that all statements made on information and belief are believed to 
be true; and further that these statements were made with the knowledge that 
willful false statements and the like so made are punishable by fine or 
imprisonment, or both, under section 1001 of Title 18 of the United States Code, 
and that such willful false statements may jeopardize the validity of the 
application, any patent issuing thereon, or any patent to which this verified 
statement is directed. 
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SYSTEM AND METHOD FOR GEOREFERENCING DIGITAL RASTER MAPS 

Cross Reference to Related Applications: 

This application contains some specification and 
figures in common with concurrently filed, copending 
applications "System and Method for Performing Flood Zone 

Certif ications" ( / filed ) and "System and 

Method for Synchronizing Raster And Vector Map Images" 

( / filed ) , which are hereby incorporated by 

reference . 

BACKGROUND OF THE INVENTION 
Technical Field: 

The present invention generally relates to geographic 
map processing and in particular to a system and method for 
associating digital raster maps with known geographic 
features. Still more particularly, the present invention 
relates to a system and method for associating specific 
points on digital raster maps with a geographic coordinate 
system. 

Description of the Related Art: 

A digital raster map is a computerized map image that 
resembles a conventional paper map in that it presents an 
image of the mapped area; but has no additional underlying 
data associated with the features of the map, A raster map 
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is typically created by scanning a conventional paper map, 
and is a grid-based map composed of pixels (or dots) of 
color or black & white. Each pixel in the grid has can be 
referenced by the pixel coordinates, and has only one value, 
5 which indicates the color of that pixel. Raster images are 

commonly referred to as "bit mapped." 

A vector map uses lines and polygons, rather than 
pixels, to describe an image. Unlike a raster map, in which 

10 the map image is simply stored as a grid of pixels, when a 

^ vector map is displayed, it is drawn from a set of 

Lq underlying data. The vector map is created through the 

fi placements of nodes on a plane and connecting those nodes 

fo with lines. Vector lines can be attributed with tables of 

ife data such as elevations, values, names or other information 

s relative to the line. Vector data can be displayed in three 

r; dimensions if the lines are attributed with z values, 

PJ modified or changed relative to user need, or layered to 

JSJ allow for turning off and on the viewing of different 

33) information. 



Because of their feature attribution properties, they 
are particularly useful for displaying geographic data. 
Vector maps are used to display boundaries or lines that 
2 5 denote the position and extent of features, such as county 

boundaries or lines denoting stream and river systems. It 
is also very easy to view or manipulate the data underlying 
a vector map, for example to view or change the elevation of 
a feature . 

30 
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Also because vector maps are commonly drawn from a 
geographic data set describing the area shown, they are very 
easily, and generally inherently, georef erenced. 
Georef erencing is the process of relating source coordinates 
to referenced geographic coordinates, which are typically in 
standard latitude/longitude. An image or a vector file is 
georef erenced to be used within a mapping/geographic 
environment. In a vector map, the data from which the map is 
drawn will typically already include a geographic coordinate 
set . 

Modern GIS systems normally make use of digital vector- 
based map information. However, a vast legacy of paper- 
based map information exists. It is very expensive and time 
consuming to convert all of the information on these paper 
maps over to a digital vector format. In many cases the 
scope and expense of such conversions render them completely 
impractical. However, even when a complete conversion to 
digital vector-based format is not possible, it is still 
possible to obtain some of the benefits of computerized map 
systems, first by converting the paper maps to digital 
raster maps (by scanning them) , and then by georef erencing 
the raster image. After georef erencing, there should be a 
clear relationship between the pixel coordinates in the 
raster map, and the geographic coordinates of the feature 
represented by that pixel. It would be desirable to provide 
a system and method for georef erencing a raster map by 
associating points on that map with corresponding points on 
a previously-georef erenced vector map, or with reference 
points designated in the raster map which have known 
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SUMMARY OF THE INVENTION 

It is therefore one object of the present invention to 
provide improved geographic map processing. 

It is another object of the present invention to 
provide an improved system and method for associating 
digital raster maps with known geographic features. 

It is yet another object of the present invention to 
provide an improved system and method for associating 
specific points on digital raster maps with a geographic 
coordinate system . 

The foregoing objects are achieved as is now described. 
The preferred embodiment provides a system and method for 
georef erencing digital raster maps by associating specific 
points on the raster map with corresponding points on a 
georef erenced vector map or another georef erenced raster 
map. According to the preferred embodiment, a raster map 
and a corresponding vector map are simultaneously displayed 
to a user. The user then locates a common geographic point 
or feature on each map, and marks each of them as a unique 
point-pair. When the user has marked at least two point- 
pairs, the system then computes a georef erencing function, 
based on the pixel-coordinates of the points marked on the 
raster map and the corresponding geographic coordinates of 
the points on the vector map. Thereafter the geographic 
coordinates of any point on the map may be easily computed. 
The preferred embodiment provides that as more point -pairs 

Fort Worth/0151MC-34350/78880 .1 



0151MC-43715 



- 6 - 



are identified, the georef erencing function is modified for 
increased accuracy. 

The above as well as additional objectives, features, 
and advantages of the present invention will become apparent 
in the following detailed written description. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The novel features believed characteristic of the 
invention are set forth in the appended claims. The 
invention itself however, as well as a preferred mode of 
use, further objects and advantages thereof, will best be 
understood by reference to the following detailed 
description of an illustrative embodiment when read in 
conjunction with the accompanying drawings, wherein: 

Figure 1 depicts a data processing system in accordance 
with a preferred embodiment of the present invention; 

Figure 2 is an image of an exemplary raster map, in 
accordance with the preferred embodiment; 

Figure 3 is an image of an exemplary vector map, 
corresponding to the raster map of Figure 2, in accordance 
with a preferred embodiment of the present invention; and 

Figure 4 is a flowchart of a process in accordance with 
a preferred embodiment of the present invention. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

With reference now to the figures, and in particular 
with reference to Figure 1, a block diagram of a data 
processing system in which a preferred embodiment of the 
present invention may be implemented is depicted. Data 
processing system 100 includes processor 102 and associated 
L2 Cache 104, which in the exemplary embodiment is connected 
in turn to a system bus 106. System memory 108 is connected 
to system bus 106 , and may be read from and written to by 
processor 102 . 

Also connected to system bus 106 is I/O bus bridge 110. 
In the exemplary embodiment, data processing system 100 
includes graphics adapter 118 connected to bus 106 , receiving 
user interface information for display 120. Peripheral 
devices such as nonvolatile storage 114, which may be a hard 
disk drive, and keyboard/pointing device 116, which may 
include a conventional mouse, a trackball, or the like, are 
connected to I/O bus 112. 

The exemplary embodiment shown in Figure 1 is provided 
solely for the purposes of explaining the invention and those 
skilled in the art will recognize that numerous variations are 
possible, both in form and function. For instance, data 
processing system 100 might also include a compact disk read- 
only memory (CD-ROM) or digital video disk (DVD) drive, a 
sound card and audio speakers, and numerous other optional 
components. All such variations are believed to be within the 
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spirit and scope of the present invention. Data processing 
system 100 is provided solely as an example for the purposes 
of explanation and is not intended to imply architectural 
limitations . 

5 

The preferred embodiment provides a system and method for 
georef erencing digital raster maps by associating specific 
points on the raster map with corresponding points on a 
georef erenced vector map or another georef erenced raster map. 
10 According to the preferred embodiment, a raster map and a 

y corresponding vector map are simultaneously displayed to a 

m user. The user then locates a common geographic point or 

feature on each map, and marks each of them as a unique point - 
m pair. When the user has marked at least two point-pairs, the 

S5 system then computes a georef erencing function, based on the 

- pixel -coordinates of the points marked on the raster map and 

Hi the corresponding geographic coordinates of the points on the 

fU vector map. Thereafter the geographic coordinates of any 

Sf point on the map may be easily computed. The preferred 

|j) embodiment provides that as more point -pairs are identified, 

the georef erencing function is modified for increased 
accuracy. 



2 5 Figure 2 is an exemplary raster map, in accordance with 

the preferred embodiment. This exemplary map shows a scanned 
image from a Federal Emergency Management Agency (FEMA) paper 
map. This raster image shows a land area with flood zone 
indications, but would, in a computer system, contain no 

3 0 underlying data regarding the area shown. 
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Figure 3 is an exemplary vector map, corresponding to the 
raster map of Figure 2, in accordance with a preferred 
embodiment of the present invention. This map shows the same 
5 area as the map in Figure 2, but is created by a computer 

system from a database describing the locations of features 
such as the streets shown. Typically, each feature shown on 
a vector map such as this will already be georef erenced, in 
that the geographic coordinates of each feature will also be 
10 recorded in the underlying data. 

Ill A digital map image is considered georef erenced if a pair 

of mathematical functions, f, and g, have been defined that 
pi can be used to convert back and forth between the coordinates 

15 of the map image (as defined by the pixels of the image) and 

s the corresponding longitude and latitude of the location of 

~r*. that point. That is, f and g do the following: 

2? 1. If (x,y) represents a location on the digital map 

S) image, then f (x,y) = (Lon, Lat) represents the longitude and 

latitude of the corresponding physical location. 

2. If (Lon,Lat) represents a physical location that 
lies within the region covered by the map, then 
25 g (Lon, Lat ) = (x, y) represents the point on the digital map image 

that corresponds to that longitude and latitude. 

Here, x and y represent the natural internal coordinate 
system of the map image. Typically, as described above, a 
3 0 digital raster map image uses the pixels of its image as a 

Fort Worth/0151MC-34350/78880.1 



0151MC-43715 



- 11 - 



natural coordinate matrix. However, in most cases, a vector- 
based map image uses longitude and latitude as its internal 
coordinate system; if so, it can be considered to be trivially 
georef erenced already. Therefore, the functions f () and g() 
above are non-trivial georef erencing functions required to 
convert back and forth between coordinate systems. 

Following is a description of a georef erencing process in 
accordance with the preferred embodiment: 

First, the raster map to be georef erenced (referred to as 
Mapl) is displayed on a computer monitor in conjunction with 
another previously georef erenced map (referred to as Map2 ) 
(step 400) . Typically, Map2 will be a vector map, but could 
also be another raster map, for which georef erencing functions 
are already known. Mapl and Map2 can be independently 
manipulated for rotation, position, and scale. The user will 
manipulate the maps until the region shown in Mapl is 
contained within the region shown by Map2 . 

Next, pairs of points representing identical geographic 
locations on Mapl and Map2 are identified, and the user uses 
a keyboard, mouse, or other pointing device to mark these 
point-pairs (step 405) . For each point-pair, one point is 
marked on Mapl and the other point is marked at the 
corresponding location on Map2 . Examples of common 

georef erencing point-pairs include street intersections, 
places where rivers cross roads, mountain peaks, building 
locations, and other readily identifiable map features. If 
any points on Mapl have points of known longitude and latitude 
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already marked on the map, then a georef erencing point -pair 
can be obtained by marking that point on Mapl and moving to 
the point with identical latitude and longitude on Map2 and 
marking it as the other half of the point -pair. 
5 Alternatively, once the known points are marked with the known 

latitude and longitude, they can be used to determine the 
georef erencing functions without requiring a corresponding 
mark on Map2 . 

10 After a point-pair is marked, the x-y point on Mapl is 

O assigned the corresponding latitude and longitude of its 

S matching point on Map2 (step 410) . 

m When two georef erencing point -pairs are marked, the 

system will compute georeference functions for Mapl based on 
s a linear transformation that allows an arbitrary rotation and 

J=* predefined scaling, between the x-y domain and the latitude- 

f|j longitude domain (step 415) . 

fj To improve accuracy and allow increased generality, more 

georef erencing point -pairs will frequently be desired. To 
facilitate the process of gathering these additional 
georef erencing point -pairs, Mapl and Map2 are now selectively 
synchronized. This means that as either map is zoomed, 

25 panned, scrolled, or otherwise caused to display a different 

region, then the other map automatically does the same. This 
greatly increases the ease with which georef erencing point- 
pairs can be identified and marked on the maps, since when the 
user of the system locates an identifiable feature on one map, 

3 0 the other map will automatically be showing very nearly the 
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same location on the other map. 

Thereafter, when the user marks a georef erencing point on 
one of the maps, then the system automatically marks the 
5 corresponding point on the other map, based on its current 

georef erencing function (step 420) . The user can then make 
corrections to the location of a point marked by the system in 
this way (step 425) . The user may enter as many point-pairs 
as are desired; each additional point-pair will increase the 
10 accuracy of the georef erencing process. 

H! When three or more georef erencing point -pairs are 

:f! determined, the system computes georeference functions for 

|0 Mapl 1 based on a completely general linear transformation 

%S between the x-y domain and the latitude-longitude domain. 

rt When four or more georef erencing point -pairs are 

Ey determined, the general linear georef erencing functions are 

2 over- determined. This means that more than the required 
O) amount of information to compute the general linear 

georef erencing functions is available, but that it is not, in 
general, completely consistent. The system use the extra 
information contained in the additional georef erencing points 
to provide validation checks to protect against the 
25 possibility that some of the data points may be inaccurate 

(step 430) . Points that deviate excessively with respect to 
a calculated standard error are presumed to be inaccurate and 
are omitted from the calculation of the georef erencing 
functions. Note that as new points are added, the system also 

3 0 rechecks points previously marked as inconsistent, to 
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determine if those points should now be considered when 
recomputing the georef erencing functions. 

The additional point-pairs allow the system to compute 
5 the general linear georef erencing functions which best fit the 

combination of all the available georef erencing point-pairs, 
so that the active georef erencing functions are revised with 
each new point-pair (step 435) . According to the preferred 
embodiment, this is done by using a "least square" parameter 
10 fitting operation. 

m The user may then proceed to enter the next point -pair 

f* (step 440) . When the user is finished, the system stores the 

fig active georef erencing functions with the raster-map (step 

Jp? 445) . At this time, the raster map is considered fully 

georef erenced. When accessed at any future time, the system 
|T§ may simply retrieve the georef erencing functions, and apply 

PJ them to find the latitude and longitude of any point on the 

;S raster map. 

£1 

The process of determining the georef erencing function 
set from a set of point-pairs is believed to be within the 
ability of one of ordinary skill in the art. The specific 
approach used by the system and method of the preferred 
25 embodiment is discussed below. 

At any given point in the georef erencing process a set of 
points p , has been specified. Let i&P denote one of the 
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points of p . Associated values x, ? y i9 Lat i!t and Lon i refer 

to the x anc * y coordinates of the point on the bitmap, and 
the latitude and longitude of the point on the earth's 
surface. We also define Ac. P to be the subset of "active" 

5 points, i.e. points which are currently being used in the 

calculation of the georef erencing functions. Assume that we 
have selected a parametric family of functions F . Functions 

in this family map f rom{x 9 y) to{Lon,Lai) . Any function f e calF 
has components fi{x 9 y) = Lon 9 and f 2 {x>y) = Lat From within this 

CI A 

3jf) family we seek the specific function, f , which (in a sense to 

yy be defined below) comes closest to making the following system 

L. of equations true: 

| f(x v y^[Lon v Lat)j for it A (1) 

O A 

Once determined, / will be the georef erencing 

15 function which is used to compute corresponding latitude and 

longitude values, i^Lon 9 Lai) for any point, (x ?> y) on the 

bitmap. There are any number of possible ways to define the 
function that "comes closest to making (1) true." We shall 
follow a "least squares" approach, also known in mathematics 

2 0 as an L 2 approach. This approach seeks to find the 

function,/, which minimizes the sum of the squared 
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differences between the actual and the predicted values of 
latitude and longitude. In other words, from among all the 

functions f e F ? f is the one which minimizes: 

SSE^[{fi^yi)~ ^^{/^{x^y)-- Lat) 2 ] (2) 

Among various alternative methods for choosing the 

A 

function / are choosing it so that it minimizes the sum of 

absolute errors (rather than squared errors, or so that it 
minimizes the largest error. Other criteria are also 
possible . 

In what follows we shall describe methods for computing 
georef erencing functions based on both a "rotational linear" 
parametric family of functions and a "general linear" family 
of functions. A knowledgeable practitioner will readily 
perceive that this approach may be easily generalized to 
other parametric families of functions. Obvious candidates 
include nonlinear parametric families of functions arising 
from the projection of the earth's surface onto a flat map. 

The General Linear Case: In the general linear case, we 
let F be the set of all possible linear transformations 
which map from (x,y) to (Lon,Lat) . Thus 
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(3) 



for some choice of parameters a ll9 a l29 <3 21? a 22 , b\> an< 3 b 2 . 

If the region covered by the map to be georef erenced is not 
too large, then this family of functions will contain a 

A 

suitable function, / whose total error is quite small. In 

the case where the map to be georef erenced covers a larger 
area than this, then the curvature of the earth must be 
taken into account and F is not a suitable family of 
functions. In such a case, nonlinear functions must be used 
as mentioned above. We shall not pursue that case further, 
since it is a straightforward extension of the procedures 
used in the linear case. 

A 

To find / we seek the parameters which minimize 

sse = X { a n x < + a uy, + b i- Lon ) 2 + { a i\ x i + a 22y< + b 2 - Lat ) 2 ( 4 ) 

The parameter values which minimize this expression are 
found by solving the following two independent systems of 
linear equations: 
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These systems can be easily solved by well known methods, 
such as Gaussian Elimination, or LU factorization. The 

solutions yield the desired values of a U9 a l29 &2\> 3 22 > 
and b 2 . 
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It should be noted that equations (5a) and (5b) do not 
have a unique solution unless three or more non-colinear 
points are contained in A . Generally speaking, then, it 
requires 3 points to choose a georef erencing function from 
the family of general linear transformations. When there 
are four points or more, it is possible to compute a 
standard deviation of error, s using the formula: 



s - 




(6) 



s is an estimator for the amount of error to be expected 
between actual and predicted latitude and longitude values. 

Note that the inverse georef erencing function, / that 
maps from (Lon, Lat) to (x,y) is readily obtained now, by 
inverting the function /. Having done this, it is possible 

to compute a similar standard deviation of error for f~ l 

which is an estimate of error as measured in the bitmap 
coordinates . 

The Rotational Linear Case: In the, so called, 
rotational linear case, we let F be the set of all possible 
linear transformations which map from {x,y) to (Lon, Lat) 
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and which also allow for a known scale change caused by the 
difference in distance per degree of latitude and per degree 
of longitude at any point not on the equator. Thus 



(7) 



10 



15 



r 

-A 



A 



y 



A, 

r 

A 



for some 



choice of parameters f$ v /? 2? /? 3? and /? 4? where 

y = cos(Lat) is computed from the estimated latitude of the map 

to be georef erenced. Clearly, this is a special case of the 
more general linear transform described above, where we take 



4? 



20 



a 22 = -fi 3 ,l\ = p x t y,b 2 = ft 2 . 



In theory, if the region covered by the map to be 
2 5 georef erenced is not too large, then this family of 

functions will contain a suitable function, / whose total 

error is quite small. In the case where the map to be 
georef erenced covers a larger area than this, then the 
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curvature of the earth must be taken into account andF is 
not a suitable family of functions. In such a case, 
nonlinear functions must be used as noted above. In 
practice, the general linear transform will yield results 
which are somewhat superior to the rotational linear 
transform due to the former's innate ability to use 
arbitrary scaling factors. These arbitrary scaling factors 
might compensate, for example, for scaling flaws in the 
digital scanner used to scan in the maps. The primary 
advantage of the rotational linear transform is that it can 
be computed with only two reference points (unlike the 
general linear transform which requires three reference 
points) . 

To find / we seek, as before, the parameters which 
minimize SSE as defined in equation (4) . 



The parameter values which minimize this expression are 
found by solving the following system of linear equations: 



n 


0 


5>. 

ieA 




0 


n 


-ly, 




I*. 


-ly, 

i<=A 


Iti + y?) 


0 


-ly, 




0 


I (*, 2 + y, 



A 
A 
A 
A 



Y J Lat 1 

1<=A 

yY J x l Lon i -Y J y,Lat l 

-rYjy\ Lon rlL x i Lat i 



ieA 



(8) 



These systems can be easily solved by well known 
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methods, such as Gaussian Elimination, or LU factorization. 

The solutions yield the desired values of j3 l9 and 

$ , which in turn yield the desired values for a u , 5 12? tf 213 
# 22? and Z? 2 . 

It should e noted that equation (8) does not have a 
unique solution unless two or more points are contained in 
A. Generally speaking, then it requires two points to 
determine a georef erencing function from the family of 
rotational linear transformations. When there are three 
points or more, it is possible to compute a standard 
deviation of error, s using the formula: 



i 



£ {a n x t + a X2 y i + b l - Lon^f + {a 2l x t + a 22 y l + b 2 - Lat^f 



n-2 



(9) 



S is an estimator for the amount of error to be expected 
between actual and predicted latitude and longitude values. 

Note that the inverse georef erencing function f that 
maps from (Lon, Lat) to (x ?< y)is readily obtained now, by 

inverting the function f . Having done this, it is possible 
to compute a similar standard deviation of error for 
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/ which is an estimate of error as measured in the bitmap 
coordinates . 

Automatic Error Detection and Handling When individual 
points are being assigned x, y 9 Lon, and Lat values, there 

is always a potential for error. To reduce the risk of 
incorrect georef erencing resulting from such errors, certain 
error handling procedures are built into the georef erencing 
process. The fundamental concept is that of detecting a 
"bad" point and then removing it from the set of active 

points, A . Note that removing a point from A will not 
delete the information associated with that point, but it 
will cause the georef erencing parameters to be completely 
uninfluenced by that point. We do not wish to remove the 
point entirely, since it may be determined at a later stage 
of the georef erencing, that the point was not really bad at 
all, and should be used in the georef erencing calculation. 
This will be clarified shortly. 

Detecting Bad Points The following steps outline the 
bad point detection process using the general linear 
transform approach to georef erencing . 

1. Begin by placing all existing points into the 
active set, A . 

2. If there are fewer than five active points then 
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you are done, Otherwise, for each of the currently active 
points in turn, move it (call it point k for the sake of 
convenience) temporarily out of the active set, and then 
calculate the resulting inverse georef erencing function 

(call it ) and its corresponding SSE k . Also, calculate 

the difference between the predicted value and the actual 



value S k 



g {k) (Lon k ,Lat k )-(x k ,y k ) 



Make a note of the 



values, <^ and 5 k I SSE k . Return point k to the active set 
(and move on to the next value of k . 

3. From among the results found in step 2. above, 
find the point, k , with the largest value ot5 k ISSE k which 

also satisfies 8 k I SSE k ) c x and S k ) c 2 , where c x and c 2 are some 

constants which are set according to the general level of 
accuracy to be expected on the particular type of map which 
is being georef erenced, the current number of active points, 
and the dots per inch of the scanned image. If there is 
such a point then mark it as bad (by removing it from the 
active set) and return to step 2 above. Otherwise you are 
done . 

There are several things to note about this procedure. 
One is that it allowing the values ofc x andc 2 to change with 

the number of active points, makes it possible for the 
georef erencing system and method to utilize points which it 
might originally determine bad or inconsistent after a large 
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enough sample of points has been gathered to make it clear 
that a lesser level of accuracy is all that can be achieved 
on this map. Another observation is that using this 
procedure it is impossible to reduce the number of active 
points down to less than four (unless you started with less 
than 4 in which case this procedure does not apply at all) . 
This scheme means that as each new point is added, all 
points determined so far are considered, even those which 
had previously been marked bad. Thus early "mis judgements" 
on the part of the system can be corrected later, in light 
of new point information. 

The same bad point detection process, can also be 
implemented using the rotational linear transform approach. 
In this case the method is capable of reducing the number of 
active points down as low as three (rather than four for the 
general linear transform approach outlined above) . This can 
be useful when dealing with small sets of active points. 

When the system is georef erencing raster map images 
that cover a large enough area so that the (nonlinear) 
curvature of the earth is a source of significant error, the 
system can calculate nonlinear georef erencing functions 
whose form corresponds to the map projection that was used 
to create Mapl . 

A specific example of the operation and application of 
the preferred georef erencing method may be shown with 
reference to the "Flood Zone Determination" business. The 
Federal Emergency management Agency (FEMA) publishes a 
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library of tens of thousands of paper maps showing various 
types of flood zones and their locations in the United 
States. A flood zone determination on a property is 
frequently done in the following way: 

1. The address of the property is examined, and the 
location of the property is determined (perhaps through the 
use of a geocoding system, or by examining an available 
street map) . 

2 . A map analyst attempts to determine which of the 
many thousands of FEMA flood maps will contain this 
property. 

3 . The map analyst goes to a map storage area and 
retrieves the desired map, often examining several maps 
before making a final selection. 

4. Having retrieved the paper map, the map analyst 
next determines where, precisely, the property is located on 
the map. 

5. Finally, the map analyst examines flood zone 
notations on the map at the property's location in order to 
determine its flood zone status. 

When performed using paper maps, the above process is 
difficult and quite time consuming. A database of scanned 
raster map images, alone, can be used to reduce the time and 
effort expended on step 3 above. However, georef erenced 
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raster map images can also be used to eliminate virtually 
all work from steps 2 and 4 above. Since those are 
generally the most time consuming steps under the current 
process, the value of georef erenced flood maps is 
5 considerable . 

Using georef erenced flood map raster images, steps 2 
and 4 above, are replaced by: 

10 2. A computer system combines the pre-designated 

Q outlines of the raster map and the georef erencing 

|K information to obtain a polygon expressed in terms of 

lU latitude and longitude that outlines the region included in 

f£ each flood map. Then the system determines which of the 

IB polygons contain the address in question, which is done 

™ using a "point- in-polygon" algorithm. At the conclusion of 

H| this process, the computer system has identified a map panel 

m (or perhaps a small number of map panels) that contains the 

address. 

5 

4. Since the latitude and longitude of the property 
are known (by virtue of a geocoding phase) , the computer 
system can use the georef erencing of the map panels to 
locate the property on each of the panels found above, thus 
25 largely eliminating any need for the map analyst to scan the 

flood map for the address location. 

Therefore, it is easily seen that by using 
georef erenced raster maps instead of paper maps, the process 
3 0 of determining flood zones is substantially automated and 
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much more efficient. 

It is important to note that while the present 
invention has been described in the context of a fully 
functional data processing system and/or network, those 
skilled in the art will appreciate that the mechanism of the 
present invention is capable of being distributed in the 
form of a computer usable medium of instructions in a 
variety of forms, and that the present invention applies 
equally regardless of the particular type of signal bearing 
medium used to actually carry out the distribution. 
Examples of computer usable mediums include: nonvolatile, 
hard- coded type mediums such as read only memories (ROMs) or 
erasable, electrically programmable read only memories 
(EEPROMs) , recordable type mediums such as floppy disks, 
hard disk drives and CD-ROMs, and transmission type mediums 
such as digital and analog communication links. 

While the invention has been particularly shown and 
described with reference to a preferred embodiment, it will 
be understood by those skilled in the art that various 
changes in form and detail may be made therein without 
departing from the spirit and scope of the invention. 
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CLAIMS : 

What is claimed is : 



1 1. A method for georef erencing a raster map image, 

2 comprising the steps of: 

3 displaying a raster map and a georef erenced map; 

4 identifying at least two geographically corresponding 

5 points on the raster map and on the georef erenced 

6 map ; 

7 associating an image coordinate of each point on the 

13 raster map with a geographic coordinate of the 
|1 corresponding point on the georef erenced map; 

M determining a functional relationship between the image 

14 coordinates and the geographic coordinates; and 
rfl thereafter, for each additional corresponding points 

13 identified on the raster map and the georef erenced 

H map, 

||| revising the functional relationship between the 

ift image coordinates and the geographic 

ff coordinates according to the additional 

18 corresponding points, and 

19 disregarding any points which are substantially 
2 0 inconsistent with the functional 

21 relationship. 

1 2. The method of claim 1, further comprising the step of: 

2 using the functional relationship to determine the 

3 geographic coordinates of features on the raster 

4 map . 



Fort Worth/0151MC-34350/78880 .1 



0151MC-43715 



- 30 - 



1 3. The method of claim 1, further comprising the step of: 

2 storing the functional relationship with the raster 

3 map . 

1 4. The method of claim 1, further comprising the step of: 

2 when the raster map is manipulated by a user, 

3 manipulating the georef erenced map accordingly. 

S3. 5. The method of claim 1, wherein the geographic 
=3 coordinates are latitude and longitude. 

s| 6. The method of claim 1, wherein the raster map and the 
4| georef erenced map are displayed on the same computer 

J 3 display. 

ml 7. The method of claim 1, wherein the corresponding points 
tfl are marked by a user after visually determining 

S geographically corresponding points. 

1 8. The method of claim 1, wherein the functional 

2 relationship is represented by a set of general linear 

3 functions . 
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4 9 . A computer system, having at least a processor 

5 connected to communicate with a readable and writeable 

6 memory , compr i s ing : 

7 means for displaying a raster map and a georef erenced 

8 map ; 

9 means for identifying at least two geographically 

10 corresponding points on the raster map and on the 

11 georef erenced map; 

12 means for associating an image coordinate of the each 

13 point on the raster map with a geographic 

fl coordinate of the corresponding point on the 

J| georef erenced map; 

M means for determining a functional relationship between 

III the image coordinates and the geographic 

#8 coordinates; and 

19 for each additional corresponding points identified on 

S) the raster map and the georef erenced map, 

M means for revising the functional relationship 

gj> between the image coordinates and the 

if geographic coordinates according to the 

24 additional corresponding points, and 

2 5 means for disregarding any points which are 

26 substantially inconsistent with the 

27 functional relationship. 

1 10. The system of claim 9, further comprising: 

2 means for using the functional relationship to 

3 determine the geographic coordinates of features 

4 on the raster map. 
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1 11. The system of claim 9, further comprising: 

2 means for storing the functional relationship with the 

3 raster map. 

1 12. The system of claim 9, further comprising: 

2 means for, when the raster map is manipulated by a 

3 user, manipulating the georef erenced map 

4 accordingly. 

Q 13. The system of claim 9, wherein the geographic 
{U coordinates are latitude and longitude. 

gt 14. The system of claim 9, wherein the raster map and the 
4i georef erenced map are displayed on the same computer 

S "S display. 

pi 15. The system of claim 9, wherein the corresponding points 

are marked by a user after visually determining 
rj| geographically corresponding points. 

1 16. The system of claim 9, wherein the functional 

2 relationship is represented by a set of general linear 

3 functions . 
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SYSTEM AND METHOD FOR GEOREFERENCING DIGITAL RASTER MAPS 
ABSTRACT OF THE DISCLOSURE 



The preferred embodiment provides a system and method 
for georef erencing digital raster maps by associating 
specific points on the raster map with corresponding points 
on a georef erenced vector map or another georef erenced 
raster map. According to the preferred embodiment, a raster 
map and a corresponding vector map are simultaneously 
displayed to a user. The user then locates a common 
geographic point or feature on each map, and marks each of 
them as a unique point-pair. When the user has marked at 
least two point -pairs, the system then computes a 
georef erencing function, based on the pixel -coordinates of 
the points marked on the raster map and the corresponding 
geographic coordinates of the points on the vector map. 
Thereafter the geographic coordinates of any point on the 
map may be easily computed. The preferred embodiment 
provides that as more point-pairs are identified, the 
georef erencing function is modified for increased accuracy. 
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